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Abstract 

Starting with a confining linear Lorentz scalar potential Vg and a Lorentz vector potential V^ 
*~~5 which is also linear but has in addition a color-Coulomb attraction piece, —Og/r, we solve the Dirac 

CN equation for the ground-state c and u quark wave functions. Then, convolving Vv with the u-quark 

'~Z^ density, we find that the Coulomb attraction mostly disappears, making an essentially linear V^ 

' ^ for the c-quark. A similar convolution using the c-quark density also leads to an essentially linear 

O 

^ Vv for the n-quark. For bound c-c charmonia, where one must solve using a reduced mass for the 

c-quarks, we also find an essentially linear V^. Thus, the relativistic quark model describes how 
j>. the charmed-meson mass spectrum avoids the need for a color-Coulomb attraction. 
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I. INTRODUCTION 

Non-relativistic potential models for charmonium states p[] have done remarkably well 
at describing these states with a simple linear confining potential (plus spin-spin and spin- 
orbit terms), absent any evidence for a short distance color-Coulomb contribution. The 
last is somewhat surprising as, early on, the high mass of the charm quark suggested that 
the color-Coulomb region might be discernible, at least in the wave functions, if not the 
eigenenergies. 

Although a Dyson- Schwinger approach would be even more appropriate |2], we report here 
on the results of attempting a self-consistent potential approach using the Dirac equation 
with both color Lorentz vector and confining color Lorentz scalar potentials designed to 
match all available data. Our principle result is that an iterative self-consistency requirement 
on these potentials leads to a color vector potential that is virtually indistinguishable from 
linear over the full range of interest. 

We take our viewpoint from the relativistic approach to the Hydrogen atom: the Dirac 
equation is used for a reduced mass electron in a potential determined by the total charge 
interior to the radial point under consideration. This approach has been studied intensively 
by, for example. Friar and Negele [3]. The closest analogy, then, is to the case of D-mesons, 
with one light quark and a charm quark in the analogous role to that of the proton in 
Hydrogen. 

As the light quark mass is negligible on the scale of interest (before we approach refine- 
ments of electromagnetic accuracy), there is no discernible reduced mass effect to consider. 
Furthermore, in Hydrogen, the charge distribution within the proton (or nucleus in more 
massive atoms) smears out the Coulomb divergence at zero separation. Here, the charm 
quark is the color Coulomb source but has no intrinsic internal structure. However, unlike 
the electromagnetic case, the virtual emission and reabsorption of gluons, produces fluctua- 
tions in the color source location that are not negligible even to leading order. Our approach 
is to take, as a flrst approximation, the charge distribution (rms size) of the D-meson as 
setting a relevant scale for these fluctuations in smearing the color Coulomb divergence. 
Although it may not be precise, given that the rms size is comparable to the inverse of the 
QCD mass scale, we consider this to be a very reasonable starting point. We use this scale to 
build a smooth truncation of the color Coulomb potential with a quadratically flat "bottom" 



and with slope equal to that of the —a^jr color Coulomb potential at the matching point. 

We have checked that such a potential produces a reasonable representation of the char- 
monium spectrum (~ ±50 MeV). We then convolve the light quark wavefunction with the 
color Coulomb potential in the manner of Friar and Negele [3] to define the potential that 
the charm quark is subject to in the presence of the light quark. Finally, to check for 
self-consistency, we again convolve this charm quark wavefunction with the color Coulomb 
potential and observe its effect on the light quark. We find a consistent, almost precisely 
linear, effective vector potential radially out all the way to the region where linear confining 
potentials are necessary for consistency with data. 

We then repeated this sequence for charmonium states where reduced mass effects are no 
longer negligible. We find a very similar result and conclude that there is indeed a single, 
consistent, approximately linear color vector potential that reasonably describes all of these 
states, despite the complications one might expect from the detailed issues described above. 
Thus, there is no remaining evidence of the color Coulomb potential despite the relatively 
large mass of the charm quark. The bottom quark may another matter. 

II. OUR ORIGINAL SET OF POTENTIALS 

After some experimentation (and with some prejudices) we came to consider the charmed 
D-mesons [c-m] and (bound) charmonia ^-c\ as resulting from scalar and vector potentials, 
Vsix) and V^{x), like those shown in Fig. 1. 

These potentials are dimensionless functions of a dimensionless radial coordinate x = nr, 
where k"^ is taken here to be 0.9 GeV/fm |1], i.e., k = 2.136 fm~^. The asymptotically linear 
slopes of Vs and V^ are the same, in accordance with the small spin-orbit splitting in the 
baryon spectrum. [5l |6] The Vg is quadratic out to x = 1.5 after which it is strictly linear: 



, x'^/(2xs) ii X < Xs, 
K(x) = <{ ^^ '' (1) 

X — Xs/2 otherwise , 

where the parameter Xs is, for us, fixed at 1.5. The flatness near x = is to preserve chiral 
symmetry at short distances. 

The Vy has, in addition to the linear confinement, a color-Coulomb contribution, 

Vy{x) = -+x-Xy, (2) 
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FIG. 1: The scalar and vector potentials Vs{x) and Vy{x) as functions of a dimensionless radial 
coordinate, x. 

as shown by the dashed part in Fig. 1. (For this figure, we used x^ = 1.0 and a = 0.4.) 
However, this is the potential seen by, say, the hght w-quark in the field of the heavy c-quark, 
which is itself moving about somewhat in the field of the w-quark. Thus it is reasonable to 
moderate the singularity at x = 0. We did this simply by altering the potential to 

(y.g {X OX f) j / Z:X ry \ X Xy^ II X "C X f) ^ 

—as/x + X — x^, otherwise , 

assuming the smoothing to be about the size of the (electric) charge radius {Rd) of the 
1^(1869) meson. Here xd = 0.16 is a reasonable guess (corresponding to Rd ~ 0.3 fm). 

It turns out that solving the radial Dirac equations [7] for a charmed quark mass of iric = 
1.550 GeV for the potentials as shown in Fig. 1 leads to a reasonably accurate description of 
the masses of the (bound) charmonia c c states - r^c, J /^J, t]'^, 'ip', Xo; and X2- [H] (For these 
charmonia calculations, the mass in the coupled differential equations must be the reduced 
mass, mc/2.) Similar calculations for c u, however, do not accurately reproduce the charmed 
D-meson masses. 



III. CONVOLVING THE COULOMB POTENTIAL FOR THE c-QUARK 



A more consistent way of moderating the Coulomb singularity is first to solve for the light 
M-quark Is ground-state wave functions [7j for the potentials Vs{x) [Eq. ([I])] and K [Eq. (|3]), 
with Xy = 1.0 and ct^ = 0.4]. Then, the vector potential that the c-quark should be subject 




FIG. 2: (Color online) The scalar potential Vs{x) and the convolved vector potential Vy{x) for the 
c-quark as functions of x. 

to is the (unrounded) V^ given by Eq. ^ modulated by the density of that w-quark. That is, 
following Friar and Negele's discussion of muonic atoms [3j, the Coulomb potential should 
be convolved with the local "charge" density defined by the w-quark Dirac wave function. 
That, together with the linear contribution, gives a new vector potential. 



Vv{x) = Qin{x)/X + Qout{x) +X - X^ 



(4) 



where 



Qm{x)/X = / x'^dx' V'ii^(x')V'„,is(x') = / x'^dx' [iplix') + Iplix')] , 

X Jo X Jo 

/•OO /"OO 

<5out(a;) = -tts / x'dx' ^l^i^{x')^u,isi^') = ~(^s / x'dx' VPl{x') + i)l{x')\ , (5) 

J X Jx 

where ipa{x) is the (real) upper component and ipb{x) the lower component. Near x = 0, 

Qin{x)/x ^ ^^^^^ / x'^dx' = ^-^^ (6 

X Jo 6 

and it never gets more negative than about —0.2 before it increases again toward zero 
like —ag/r. As for Qout{x), since the Is upper component radial wave function is well- 
approximated as a Gaussian JUIH], its integral gives, approximately, — a^ times a (narrower) 
Gaussian. 

A plot of Vy{x) calculated from the integrals of Eq. (pi) is given in Fig. 2. Despite its 
appearance, it is not strictly a straight line - there is some small curvature in the plot below 
a; = 1. Nonetheless, we consider the high accuracy of a linear approximation to be rather 
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i/A (I5) for original V (x) 
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FIG. 3: Comparing the c-quark radial wave functions for the original and convoluted potentials. 

surprising, as we were expecting only a minor change in the effective value of xd- The 
dashed line in Fig. 2 is a linear fit to K)(a^), aox + ai, with slope oq = 1.139 and displacement 
ai = -1.498. 



IV. HOW Vc CHANGES WITH THE NEW POTENTIAL 

If one solves for the Is state of the c-quark for rric = 1.550 GeV and the original potentials 
of Eq. ([I| and Eq. ([s]), one finds the c-quark eigenenergy to be Ec{ls) = 1.339 GeV. That 
is, the energy of the c-quark for these potentials is some 200 MeV less than its mass. The 
upper and lower Is radial wave functions are displayed in Fig. 3 as dashed curves. The rise 
of ipa at the origin is reminiscent of the non-relativistic ground state wave function for a 
pure Coulomb potential, which is a wave function that is a decaying exponential like e~"^. 

However, solving for Ec{ls) with Vs{x) and the convolved Vv{x), we find Ec{ls) = 1.534 
GeV, now considerably higher in energy because of the missing Coulomb well, although still 
less than rric- The upper and lower component wave functions, shown as the solid curves in 
Fig. 3, are somewhat broader than those found for the original potentials. (In both cases, 
however, the c-quark wave functions are not as broad as those for the w-quark.) 

Thus, while the mass of a cu D-meson, which depends strongly on Ec{ls), is not very 
sensitive to the difference between the original and convolved potentials, the wave functions 
are quite different. We therefore expect that quantities such as transition strengths, which 
are more dependent on the details of the wave functions, will be more potential-dependent. 




FIG. 4: (Color online) The scalar potential Vs{x) and the convolved vector potential Vy{x) for the 
li-quark as functions of x. 

as usual. 

V. CONVOLVING FOR THE n-QUARK POTENTIAL 

The singular Coulomb piece of the vector potential seen by the tt-quark is also smeared 
by the motion of the somewhat more confined, slower moving c-quark. With formulae like 
Eqs. Q and (jSj), but with the density provided by |-?/'c(ls)P, we again found this second 
convolution is also very close to linear, as shown in Fig. 4. In making this plot we used 
the ipd^s) found from using Vy. Again, the dashed line is a linear fit to Vy{x) with slope 
ao = 1.218 and displacement ai = —1.673. 

The K) is slightly more curved and a bit deeper than the V^ shown in Fig. 2, but is 
essentially the same nearly-linear potential as that which affects the c-quark. Although 
the linear fit parameters are slightly different, one could reasonably assume the same linear 
potential for both quarks. 

VI. CONVOLVING FOR CHARMONIUM 



Similarly, we investigated the smearing of the Coulomb potential for the charmonia [c-c] 
states. In this case we must solve for the reduced mass, which is mc/2 = 0.775 GeV, but 
otherwise the calculation proceeds much as above. First we find the wave functions for the 




FIG. 5: The scalar potential Vs{x) and the convolved vector potential K)(x) for charmonia states 
as functions of x. 

c-quark for the original potentials, Vs{x) [Eq. ([I])] and K [Eq. (|3])]. Then, convolving the 
singular Coulomb potential with |'?/'c,reduccdP as in Eqs. (0) and rtsj we obtain the plot of, 
now, Vy shown in Fig. 5. The dashed line is a linear fit to Vy{x) with slope Oq = 1.225 and 
displacement ai = —1.688. Note that the linear fit parameters are extremely close to the 
values that we found for the (effective) potential for the u-quark. 



VII. DISCUSSION 

These results are consistent with the good spectral results found in non- and relativistic 
models for the charmonium spectrum using linear confining potentials. It is unclear whether 
the slight differences in the effective linear potentials merit the complications of a relativistic 
approach to the calculation of spectra until high accuracies become necessary. Our convolu- 
tion approach suggests that the good spectral results with a non-relativistic, linear potential 
are due to the still relatively light mass of the charm quark. This in turn invites the question 
as to whether this status can still hold for the bottom quark. Absent detailed comparisons 
with transition rates for charmonium and D-meson decays, which one expects to be more 
sensitive to wave function details than are spectra, the color Coulomb contribution to the 
effective potential for quark binding remains undetermined from charmonium data alone. 
We intend to turn next to bottom quark states to examine whether the spectra there can 
provide a definitive determination of this issue. 
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